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ABSTRACT 

It has been shown that resonant absorption is a robust physical mechanism to explain the observed 
damping of magnetohydrodynamic (MHD) kink waves in the solar atmosphere due to naturally oc- 
curring plasma inhomogeneity in the direction transverse to the direction of the magnetic field. Theo- 
retical studies of this damping mechanism were greatly inspired by the first observations of post-fiare 
standing kink modes in coronal loops using the Transition Region And Goronal Explorer (TRACE). 
More recently, these studies have been extended to explain the attenuation of propagating coronal 
kink waves observed by the Goronal Multi-Channel Polarimeter (CoMP). In the present study, for the 
first time we investigate the properties of propagating kink waves in solar waveguides including the 
effects of both longitudinal and transverse plasma inhomogeneity. Importantly, it is found that the 
wavelength is only dependent on the longitudinal stratification and the amplitude is simply a product 
of the two effects. In light of these results the advancement of solar atmospheric magnetoseismology 
by exploiting high spatial/temporal resolution observations of propagating kink waves in magnetic 
waveguides to determine the length scales of the plasma inhomogeneity along and transverse to the 
direction of the magnetic field is discussed. 

Subject headings: Sun: oscillations — Sun: corona — Sun: atmosphere — Magnetohydrodynamics 
(MHD) — Waves 



1. INTRODUCTION 

Resonant absorption, caused by plasma inhomogene- 
ity in the direction transverse to the magnetic field, has 
proved to be the most likely candidate for explaining the 
observed attenuation of magnetohydrod ynamic (MHD) 
kink waves in the sola r atmosphere (see ICoossens et all 
[20061: IGoossensI [20081 for a review about this damp- 
ing mechanism). E.g., resonant absorption is a feasi- 
ble explanation for the damping of kink MHD waves 
in both coronal loops (Goosscns ct al. 2002) and in the 
fine structure of solar promin ences (see, e.g., I Oliver! 
120091 : [Arregui fc Ballestedl2010L and references therein). 
Most of theoretical studies of resonant absorption as- 
sume that the inhomogeneity is in the radial direction 
exclus ively, and only a few works (e.g., Andric s et al.l 
2005'; 'Arregui ct al. 2005; iDvmova fc Rudermann2006l : 
Andries et al.. , 2009a: Solcr e t al.ll2010t ) have studied the 
properties of resonantly damped standing kink waves 
when plasma inhomogeneity is also present in the di- 
rection of the magnetic field. These studies have been 
restricted to standing waves, so the amplitude of the 
standing wave is damped in time as a result of the reso- 
nance. Nevertheless, there are observational indications 
of spatial damping of propagating waves along differ- 
ent waveguides in the solar at mosphere. Som e clear 
exa mples have been reported bv [Tomczvk et al.l (|2007D 
and iTomczvk fc Mclntoshl ()2009D in coronal loops us- 
ing Coronal Multi-Channel Polarimeter (CoMP) data. 
Other observational evidence of propagating kink waves 
in wave guides of the solar atmosphere are, e.g., 

lroberto.soler@wis.kuleuven.bel 



De Pontieu et a l.' (2007^ : lHe et al. I (|2009al [bl) in spicules, 
Okamoto et al.l (.2007,') in the fine structure of promi- 



nences. and iLin et all ()2007l I2Q09D in filament threads 
The theoretical modeling of the spatial damping of trav- 
eling kink wav es due to reson a nt ab sorption has been 
carried out by ITerradas et al.' (2010) and Vcrth ct alj 
(jlOlO) analytically, and by Pascoc ct al. (2010, 2QX^ 
numerically. All these works assumed a homogeneous 
densit y in the longitudinal direction. ITerradas et al.l 
(|2010f ) showed that the damping length by resonant 
absorption is inversely proportional to the wave fre- 
quency. It was predicted that high-frequency waves be- 
come da mped in lengt h scales smaller than low-frequency 
waves. iVerth etaH (pOlOi) showed that this theory 
is consistent w i th th e Co MP observations repo r ted b y 
ITomczvk etall (|2007[ ) and ITomczvk fc Mclntoshl (gOQl). 
This means that in solar waveguides resonant absorption 
acts as a natural low-pass filter. 

Some rec ent investigations have extended the analyti- 
cal work of ITerradas et al.l (|2010[ ) b y adding eff e cts not 
considered in their original paper. iSoler et all (|2011aD 
took into account the effect of partial ionization in the 
single-fluid approximation, so kink waves are damped 
by both resonant ab sorption and ion-neu tral collisions in 
their configuration. ISoler et"al1 ()2011al ) concluded that 
in thin tubes the effect of re sonant absorpt i on do mi- 
nates, recovering the result of ITerradas et al.l (|2010( ) in 
the fully ionized case that the damping length is in- 
versely proportional to_the frequency. Therefore, the re- 
sult of [Soler et aI7 (2011a') allows the theory developed 
bv ITerradas et al. (2010) to be applied to kink waves in 
partially ionized plasmas as, e.g., chromospheric spicules 
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and solar prominences. Alternately, iSoler et alj ()2011bO 
have shown that in the presence of flow the damping 
length remains inversely proportional to the frequency, 
but the factor of proportionality is different for forward 
and backward propagating waves to the flow direction 
and depends on the characteristics of the flow. However, 
none of these previous works investigated the effect of 
longitudinal stratification. 

The analysis of kink MHD waves in longitudinally ho- 
mogeneous models make some problems more tractable 
from an analytical point of view. Realistic models of 
solar waveguides should include longitudinal stratifica- 
tion. E.g., coronal loops can extend out into the atmo- 
sphere up to heights of the order of several density scale 
heights. This means that the differences in the density 
at their footpoints and apex can be significant. There 
is evidence of such stratification in coro nal loops both 
from emission measure analysis (e.g.. lAschwanden et aTl 
19991) and coronal magnetoseismologv fsee lAndries et all 
2009bl for review of this topic). Regarding spicules, there 
have been spectroscopic line intensity studies to show the 
plasma density is longitudinally strongly stratified (e.g., 
[Beckers 1968; Makita 2003) and more recently by mag- 
netoseismology (IVerth et aLll2011f). It is also clear from 
Ha observations fe.g.. iLin et al.l 20081 ) that prominence 
threads too exhibit longitudinal density inhomogeneity. 
For these reasons, the main aim of the present paper is to 
investigate the effect of longitudinal stratification on the 
properties of resonantly damped propagating kink waves. 
This is an interesting problem from a theoretical, as well 
as observational point of view. The reason is that a radial 
variation of density causes damping due to resonant ab- 
sorption while a longitudinal variation might cause an 
increase of the amplitude. An increase of the ampli- 
tude due to stratification along the m agnetic field direc- 
tion has been obtained by, e.g.. iDe M oortcl ct al. (199^ 
for ph ase mixed Alfven waves and Ibe Moortel fc Hoodl 
(|2004f ) for propagating slow modes. 

In this paper, we investigate the effect of longitudinal 
density variation on propagating kink waves by using 
both analytical and numerical approaches. Analytical 
theory is used to the largest possible extent. To do so, 
we apply standard approximations usually adopted in the 
preceding literature for the investigation of kink waves in 
magnetic tubes. The approximations used in the analyt- 
ical part of this paper are briefiy discussed and justified 
in the next two paragraphs. In addition to the analyti- 
cal theory, the problem is solved numerically beyond the 
limitations of the analytical approximations. The com- 
parison between the analytical results and those from 
the full numerical simulations will enable us to check the 
validity of the analytical theory. 

Our waveguide model is a straight cylindrical mag- 
netic flux tube, inhomogeneous in both the radial and 
longitudinal directions, and embedded in a magnetized 
environment. To investigate propagating kink waves an- 
alytically we adopt the following approximations. We 
use the thin tube (TT) approximation, so that the wave- 
length, A, is much larger than the radius of the tube, R. 
This condition is easily fulfilled for propagating waves 
observed in the solar atmosphere (see the various refer- 
ences given above) . The variation of density in the radial 
direction is confined to a thin transitional layer, and we 
adopt the thin boundary (TB) approximation (see de- 



tails in, e.g.. iHoUwee & Yanslll988l 


Sakurai et al.l 119911: 


iGoossens et al.lll992l: lGoossensll2008 


). The TB approxi- 



mation has been proved to be very accurate even when 
radi al inhomogeneity is not restrict e d to a thin layer (see , 
e.g.. iVan D oorsselaere et al. 2004; Andries e t al.l 120051 : 
lArregui et al.„ .2005. ) . Regarding inhomogeneity in the 
longitudinal direction, we assume that the kink mode 
wavelength is shorter than the longitudinal inhomogene- 
ity length scale, so that the WKB approximation can be 
applied. In the numerical part of this paper, the prob- 
lem is solved in the general case where large variations 
of density in both the radial and longitudinal directions 
are considered and for arbitrary wavelength. 

Curvature is neglected in the model. The effect of 
curvat ure was studied by, e.g., i Van Doorsselaere et al.l 
(l2005f) and Tcrradas ct al. (200^), and it has been re- 
viewcd by Van Doorsselaer e et all (j2009D . Analytically, 
IVan Doors selaere ct al. (:2005[) showed that curvature has 
no first-order effect on the frequency and the damping of 
kink modes in curved cylindrical models. The result that 
curvature has only a minor effect on kink MHD waves was 
also numerically conflrmed by IVan Door sselaer e et al.l 
(^005) and Tcrradas ct al. (2dO(J^ Wave damping due 
to leakage is not considered. Leaky waves, i.e., waves 
damped by MHD radiation, wer e first studied in mag- 
netic flux tubes by iSpruiti ()1982D , whose conclusion was 
that leakage is only important for wavelengths of the 
order of the tube radius or smaller. For thin tubes 
leaka ge is unimportant. T his r esult was confirmed by 
IGoossens fc Hollwed (|1993D and IGoossens et al.i (.20091 
Finally, the /3 = approximation is also adopted, where 
/3 is the ratio of gas pressure to magnetic pressure. 
The assumption of /3 = implies that waves do not 
have motions along the equilibrium magnetic field di- 
rection. It has been ex plained by, e.g.. Spruit (1982) and 
IGoossens et al.l (j2009( ) that the longitudinal component 
of the velocity is proportional to the transverse compo- 
nent, with a factor of proportionality that depends of (i 
and {R/Xf' . In the low-/3 case and for thin tubes the lon- 
gitudinal component of the velocity is much smaller than 
the transverse component. In magnetic structures of the 
solar atmosphere /3 ^ 1. This means that the longitudi- 
nal component of the velocity can be neglected and the 
/3 = approximation can be safely adopted. For the ef- 
fect of plasma /3 o n the res onant damping of kin k waves, 
ISoler et al.l (|2009( ) (see also IGoossens et al.l[2009 ) showed 
that the contribution of the slow continuum damping, 
present when /3 ^ 0, is negligible compared to the Alfven 
continuum damping. 

This paper is organized as follows. The description of 
the model configuration and basic equations are given 
in Section [2] Then, we obtain general expressions for 
the amplitude and wavelength of propagating kink waves 
in longitudinally and transversely inhomogeneous mag- 
netic flux tubes in Section [31 Later, we use these gen- 
eral expressions to study resonantly damped propagat- 
ing kink waves in stratified solar waveguides and com- 
pare the analytical predictions with full numerical time- 
dependent simulations in Section [4l The implications 
for the method of solar atmospheric magnetoseismology 
are given in Section [SJ Finally, Section [51 contains the 
discussion of the results and our conclusions. 
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2. MODEL AND BASIC EQUATIONS 

In the (3 — approximation, the basic equations for 
the discussion of linear ideal MHD waves are, 



dv 1 
ot /i 
9b 

at 



(V X b) X B, 
: V X (v X B) , 



(1) 
(2) 



where p is the plasma density, B is the equilibrium mag- 
netic field, V is the velocity perturbation, b is the mag- 
netic field perturbation, and p is the magnetic permit- 
tivity. We use linear theory in our investigation be- 
cause the Alfven velocity in the solar corona is of the 
order of 1000 km s~^, while the peak-to-peak veloc- 
ity amplitude of the propagating waves observed using 
CoMP is only about 1 km s ~^ (see lTomcz vk ct al. H2007t 
pTomc zyk fc Mclntoshll2009f l. We study MHD waves in 
a static background, so the equilibrium fiow is absent 
from our discussion. The effect of fiow on resonantly 
damp ed propagating kink waves has been recently stud- 
ied bv lSoler et ai] (12011b '). 

The equilibrium configuration is a straight cylindri- 
cal magnetic fiux tube of average radius R embedded 
in a magnetized plasma environment. For convenience, 
we use cylindrical coordinates, namely r, and z for 
the radial, azimuthal, and longitudinal coordinates, re- 
spectively. The axis of the cylinder is set along the 
z-direction. In the following expressions, subscripts i 
and e refer to the internal and external plasmas, respec- 
tively. We denote by pi{z) and Pe{z) the internal and 
external densities, respectively. Both of these quantities 
are functions of z. There is a nonuniform transitional 
layer in the transverse direction that continuously con- 
nects the internal density to the external density. The 
thickness of the layer is I, and extends in the interval 
R — 1/2 < r < R + 1/2. The equilibrium magnetic field 
is straight and homogeneous, B = Be^, with B con- 
stant. As if is an ignorable coordinate in our model and 
we consider waves propagating along the tube with a 
fixed frequency, we write all perturbations proportional 
to exp (imif — iujt), where m is the azimuthal wavenum- 
ber and ui is the wave (angular) frequency. The frequency 
/ is related to a; by a; = 2nf. Kink modes, i.e., the 
only wave modes that can displace the magnetic cylinder 
axis and so pro duce tran sverse motions of the wh ole flux 
tube (see, e.g., lEdwin fc Roberts 1983; Goosse ns et al.l 
I2009f ). are described by m = 1. Due to the presence 
of a transverse inhomogeneous transitional layer, wave 
modes with m ^ are spatially damped by resonant 
absorption. In the position within the transverse inho- 
mogeneous transitional layer where the wave frequency 
coincides with the local Alfven frequency, i.e., the reso- 
nance position, the wave character changes and becomes 
more Alfvenic as the wave propagates. As a result of this 
process, transverse motions of the flux tube are damped 
while azimuthal motions within the transitional layer are 
amplifled. 

We use the TT approximation, i.e., A ^ i?, where 
A is the w avelength along the tube, and app ly the for- 
malism of iDvmova fc RudermanI (|2005l [200l . The ef- 
fect of the resonant damping is included here by intro- 
ducing jump conditions for perturbations at the reso- 
nance position along with the TB approximation (see de- 



tails in, e.g.. iHpUweg fc Yanel ll988t iSakurai et all 119911: 
iGoossens et al.lll992f) . The combination of both TT and 
TB approximations is usua lly called the TT TB approach 
(see the recent review by IGoossensl[200l. Then, fol- 
lowing the scaling arguments of iD vmova fc Rudernianl 
(200 5, 2006), the total pressure perturbation P = Bbz/li, 
with bz the longitudinal component of the magnetic field 
perturbation, evaluated at both sides of the transverse 
transitional layer can be cast as 

Pi«(-j A{z), at r«i?--, (3) 

/ D\ ™ r 

Pc«f-j A,{z), at r«i?+-, (4) 

with Ai(z) and Ac{z) arbitrary functions of z. In the TB 
approximation, the condition of P continuity acr oss the 
resonant la yer (e.g., Sakurai ct al. 1991; Goossens et al.l 
[1991 [T995[) implies Pi = P^ at r = R, which means that 

A{z) = Ae(z) = Aiz). 

On the other hand, from Equations ([T]) and ([2]) we 
obtain the following relation 



1 d'^P 



(5) 



dt^ "^'"'922 p{z)drdt' 

where Vr is the radial velocity perturbation, and v\{z) = 
is the Alfven velocity squared. After using Equa- 



fj.p{z) 



tions ([3]) and ^ in Equation ([5]) we arrive at the two 
following expressions. 



vliiz) + 



-^^^—A{z), 



Pi{z) R 



at 



R 



Pc{z) R 



A{z), at rKiR + 



where Vri and the radial velocity perturbations 

at r « i? — 1/2 and r Ki R + 1/2, respectively. Now, 
we take into account that in the TB approximation the 
jump of Vr across the resonance layer, namely [vr], is 
assumed to be the same as the jump of Vr across the 
whole transverse inhomogeneous layer, i.e., Vre — Vri = 
[vr]. For a straight and constant magnetic field, \vr] is 
given , by (see, e.g., Sakurai et al. 1991; Goossens ct alj 
1991 [19951: iTirrv fc Goosse ns 1996: Andries et al...2005r 



Dvmova fc Rudermanll2006l ). 



UJ \drp\ 



at 



(8) 



where \drp\ is the radial derivative of the transverse den- 
sity profile at the resonance position, which is assumed 
to take place at r « i?. This is a reasonable assumption 
in the TTTB approximation. Subsequently, we combine 
Equations ([6]) and ([7]) to eliminate function A{z) and use 
the jump condition (Equation (|8])) to arrive at a single 
equation for Vri and [vr\- Forthwith, we denote Vri as Vr 
to simplify the notation and the resultant expression is 



UJ 



d'^Vr 



vl.{z) 



dz^ v\^{z) 



where 



vl{z) 



2B' 



fJ-iPiiz) + Pe{z))'' 



(9) 



(10) 
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is the squared of the kink velocity, i.e, the phase velocity 
of the propagating kink wave. For future use, it is in- 
structive to note that the kink velocity vy^ is a function 
of z. Finally, we use Equation (j6|) along with the result 
that P is constant across the resonant layer to relate [vr] 
with Vr, obtaining 

^^"-^ R\drp\ \ ^2 0,2 +^rj , (iij 

which we put in Equation ([9]) to get an equation involving 
Vr only, namely 

with 

C4 = -i^ 



C3 = - 



mpi(z) 9/1 



C, = 1 - ^ 



vr m 
'2^ 



Piiz) 



^'Ai(^) 
1 



92^ \\drP 



Pi{z)+Pe{z) 
\drp\ 



. m d 

R oz 



\drp\ 



2 i? 



dz-^ 



Pijz) 

\drp\ 



Pi{z) 



(13) 



Equation (jl2p is the main equation of this investigation. 
Note that the imaginary parts of the coefficient functions 
C0-C4 are all proportional to m. 

2.1. Case Without Longitudinal Stratification 

As a limiting case of Equation (fT2|) . we study wave 
propagation in a magnetic tube which is uniform in the 
longitudinal di r ection and try to recover the results of 
iTerradas et"al] (|20ia hereafter TGV). In order to do so, 
we take pi and pe constants and neglect their derivatives 
in the longitudinal direction. Since the background is in- 
dependent of z we can Fourier-analyze in the z-direction 
and write Vr proportional to exp (ikzz), with the lon- 
gitudinal wavenumber. Equation (1121) becomes 



C4kt - C2kl 



Co = 0. 



(14) 



Since we are investigating spatial damping, we write 
kz = fczR + ikzi, approximate fc^j^ « w^/w^, and assume 
fczi <^ fczR, (weak damping), so we neglect terms with fc^j. 
Finally the following relation is obtained. 



kzl _ TT TO 1 (pi - po)^ 
kzR 8 R \drp\ Pi + Pc 



(15) 



which coincides with Equation (10) of TGV. Therefore, 
we recover the results of TGV for propagation and damp- 
ing of kink waves in longitudinally homogeneous tubes. 
The reader is referred to TGV for a complete analysis of 
this case. It is worth noting that Equat i on ([151) also 
agrees with Equation (20) of iSoler et al.l (|2011al ). ob- 
tained for the spatial damping of propagating kink waves 
in partially ionized threads of solar prominences, if the ef- 
fect of damping by ion-neutral collisions is removed from 
their expression. 



3. GENERAL FORMULATION 

3.1. No Resonant Damping 

In the present and next Sections, we investigate the 
effects of longitudinal density stratification. First, we 
study the case when resonant damping is absent to assess 
the effect of density stratification on the amplitude of 
the propagating kink wave. Damping will be considered 
later. Then, Equation (|12p reduces to 



d^Vr 



dz^ vl{z) 



0. 



(16) 



Note that Equation (fT6| is independent of the azimuthal 
wavenumber, to. In the absence of damping, the waves 
are degenerate with respect to m as long as to ^ 0. 

To find a solution to Equation (|T6l) . we use the 
Went zel-Kramers-Brillouin (W KB) approximation (see, 
e.g., iBender fc Orszagl 119781 for details about the 
method). We define the dimensionless longitudinal co- 
ordinate as, 

C=|, (17) 

where A is the longitudinal inhomogeneity scale height, 
and express the solution of Equation ([T5| in the following 
form, 

«r.(C)«Q(C)cxp[*A/C(C)], (18) 

with Q (C) and K. (C) functions to be determined. The 
validity of the WKB approximation, and so of Equa- 
tion (fT8|) . is restricted to large values of A, i.e., weak 
inhomogeneity, such as A/ A ^ 1. We combine Equa- 
tions ([T5| and (fTS|) . and separate the different terms ac- 
cording to their order with respect to A. As A is large, 
the dominant terms are those with the highest order in 
A. From the leading term, i.e., that with ©(A^), we 
obtain the following expression. 



where fcz(C) has been defined as 

d/c(C) 



kziC) 



dC 



(19) 



(20) 



We see that kziC) plays the role of the longitudinal 
wavenumber, which in our case is a function of C through 
the dependence of the kink velocity on ( (Equation ([TOl) ). 
Therefore, the wavelength, A(C), is 



A(C) = T 



27r 



kzRiO 



—MO 



TWk(C), 



(21) 



with T the wave period. We obtain that A(C) changes in 
C due to stratification. Therefore, function /C(C) is 



/c(C) 



c 



1 



(22) 



/o ^'k(C) 

where we have assumed /C(0) — 0. 

On the other hand, from the contribution of O (A^) in 
the WKB expansion of Equation (fTHt , we get the follow- 
ing expression 



dQiC) 1 dkziO 
dC 2kz{C) dC 



QiO = 0, 



(23) 
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which can be easily integrated to 

' \ 1/2 



g(c) = Q(o) 



^^k(0) 



(24) 



Equation (|24|) teUs us how the ainphtude changes with 
height due to density stratification. Density stratification 
modifies both the ampUtude and the wavelength of the 
kink wave as it propagates along the magnetic tube. If 
density decreases with height, and hence the kink veloc- 
ity increases with height, Equations p4)) and (|21l) show 
that both the wavelength and the amplitude increase. 

The complete expression for in the WKB approxi- 
mation using the original variable z is 



Vr{z) K. Vr{Q) ( I CXp [ lUJ 



1 



Wk(5) 



dz 



(25) 



3.2. Effect of Resonant Damping 

Hereafter, we include the effect of damping by resonant 
absorption and we try to solve the full Equation ([T^ . As 
a completely general solution to Equation (IT2|) is very dif- 
ficult to obtain, we restrict ourselves to the case in which 
the internal and external densities are proportional. We 
define the density contrast, x, as 



X 



(26) 



We assume that parameter x is a constant independent 
of z. In addition, we write \drp\ as 



\drp\ = T 



TT^ Pijz) - pc{z) 

4 / 



J" 



1 pi{z) 



I 



(27) 



with J' a factor that depends on the form of the trans- 
verse d ensity profile. For exa mple, F = A/tt^ for a linear 
profile (iGoossens et al.ll2002f) and = 2 /7r for a sinu- 
soidal profile (jRuderman fc Robertil2002l ). Then, coeffi- 
cients C0-C4 in Equations (fTS]) simplify to 



C4 = - 

C2 = - 
Ci=( 



2 / m 
-i — 



X 



1 



Aim X 1 dv\^{z) 
TT i? X — 1 dz 



.2 I m 



X I d^vU^) X + 1 
X — 1 dz^ X ~ 1 



2 I m 1 



v^{z) TT R F X - ^ vXiiz) 



(28) 



As before, we use the WKB approximation to solve 
Equation with the coefficients given in Equa- 

tions (|28p . Then, we re-introduce the dimensionless lon- 
gitudinal coordinate C, = z/K (Equation pT|) ) and write 
Wr(C) in the form given in Equation (IT51) . Following the 
same process as in the case without resonant damping, 
we separate the different terms of Equation (IT^ accord- 
ing to their order with respect to A. From the leading 
contribution, i.e., that with O (A**), we obtain the fol- 
lowing expression. 



.Trm ^ Pi(C)Pc(C) 
2 RB^ \drp\ 



= 0. (29) 



where kz{C) is defined in Equation pO)). For conve- 
nience, we have used Equation ([27| to introduce \drp\ 
again, in order to illustrate that Equation (|29|) is for- 
mally equivalent to Equation (6) of TGV if the longi- 
tudinal dependence of the density is added to their ex- 
pression. We follow the analytical process described in 
TGV and write k^iC) = fczR(C) + ik^iiC). The effect of 
resonant absorption is to introduce an imaginary part 
of fcz(C)j while we can reasonably assume that fczr{(C) 
is the same as in the case without resonant damping, 
i.e., fczR(C) ~ w/uk(C) according to Equation (IT^ . Next, 
we assume weak damping, i.e., kzi{C) <S fczR(C)j so we 
neglect terms with fc^ilC) i^i Equation (|29p . Finally, a 
relation between fczR(C) and kzi{C) is derived, namely 



fczl(C) ^™ 1 (Pi(0-Pe(C))' 



fczR(C) 8R\drp\ Pi(C)+Po(C) 



(30) 



which coincides with our Equation (jl5l) and with Equa- 
tion (10) of TGV if the dependence of the densities in 
the longitudinal direction is explicitly included. Taking 
into account the definitions of Equations (pS)) and (P71) . 
we rewrite Equation ([50]) as 



fczi(C) ^ I ml x-l 



(31) 



and we clearly see that the right-hand side of Equa- 
tion ([3T|) is independent of C- We use Equation ([2T|) 
and define the damping length by resonant absorption as 
Ly){() = l/fczi(C), so we can obtain from Equation pip 
that 

LuiO^^Rx+l , . 

A(C) mlx-l ^ ' 

Again, note that the right-hand side of Equation (|32p 
is independent of C. Finally, we consider the expression 
for A(C) of Equation (PT|) to give the expression for the 
damping length as 

iD(C) = 2^-y^it,k(C), (33) 

TO i X ^ 1 

which is equivalent to Equation (22) of TGV. The damp- 
ing length is inversely proportional to both uj and m. By 
using Equations ([21]) and ([33]) . we compute the real and 
imaginary parts of /C(C) = ^r(C) + *^i(C)i namely 



/Cr(C): 



« 2n rC 



A(0 



1 



^'k(0 



1 m I X ^ 

^2^yRx + i^ 



(34) 



(35) 



where again we assumed /Cr(0) = /Ci(0) = 0. 

On the other hand, we get an equation for Q{C) from 
the contribution of O (A^) in the WKB expansion of 
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Equation (|12l) . namely 



Equations (IMj) and (pHl) . the square of the kink velocity 



gQ(C) 1 dkMO 



QiO = 0, 



(36) 



which coincides with Equation (|23p . Therefore, the solu- 
tion to Equation p6l) is given by Equation (I24[) . mean- 
ing that the change of the amplitude in C due to density 
stratification is the same as in the case without resonant 
damping. 

Finally, using Equations ([M]), and in Equa- 
tion (|18p . and returning to the original variable z, we 
obtain 



Vr{z) ~ -4(z) exp ILJ 



1 



l'k(2) 



-dz 



(37) 



with A{z) the amplitude as a function of z given by 



Aiz) = vr{0) 



X exp 



f ^\ 1/2 

MO). 
1 mix 
2^ y^X+ l' 



1 



-dz 



(38) 



Equation psp shows that there are two effects that de- 
termine the amplitude of the propagating kink wave, i.e, 
longitudinal density stratification and resonant absorp- 
tion. While the effect of longitudinal stratification de- 
pends on whether the density increases or decrease with 
z, the role of resonant absorption is always to damp the 
wave and so to decrease its amplitude. It is clear from 
Equation (1381) that resonant absorption acts as a natural 
filter for waves with high frequencies and large m. 

4. APPLICATION TO PROPAGATING KINK 
WAVES IN STRATIFIED CORONAL LOOPS 

4.1. Analytical Expressions 

Here, we apply the general theory described in the last 
Section to the case of driven kink waves in a stratified 
fiux tube. For a particular application we model a coro- 
nal loop being driven at one footpoint. Assuming the 
longitudinal stratification of the coronal loop is symmet- 
ric about its apex, we focus our analysis on the loop half 
that contains the footpoint driver. We take the plasma 
to be more dense at the loop's footpoints to mimic a 
gravitationally stratified solar atmosphere and assume 
an exponential z dependence in density, namely 



Pi(z)=Pioexp , 



Po(z) = poocxp (--j , 



(39) 
(40) 



where pio and peo a-re the internal and external densities, 
respectively, at the footpoints of loop, i.e., at z = 0, and 
A is the scale height defined as 



A = 



L/2 



In (pio/papox) 



with L the total length of the loop and Papcx the internal 
density at the apex, i.e., at z = L/2. As we assume the 
same stratification within the tube and in the corona, 
X = Pi{z)/pe{z) = Pio/peo is a constant. According to 



(42) 



where w^q = 2B'^ / [/i (pio + Pco)] is the square of the 
kink velocity at z = 0. As in our model the mag- 
netic field is straight and constant, the density scale 
height and the kink velocity scale height are propor- 
tional. CoMP obser vations of coronal propagating waves 
(jTomczvk fc Mcinto sh 2009) show that the kink velocity 
is relatively constant with height so that the kink veloc- 
ity scale height is much longer that the typical coronal 
scale height. The effect of flux tube expansion might be 
responsible for keeping the kink velocity relatively con- 
stant. The effect of flux tube expansion is discussed in 
Section m 

4.1.1. Case Without Resonant Damping 

First, we analyze the case without resonant damp- 
ing. The kink mode propagation is governed by Equa- 
tion ([T6|) . The expression for w^(z) given in Equation (|42p 
enable us to rewrite Equation (jl6|) as a Bessel Equation 
of order 0. Then, it is possible to obtain an exact solution 
of Equation p6)) as 



Vr{z)^AiHl,'^ 



with H^^"^ and Hq^' the Hankcl functions of the first and 
second kind of order 0, and Ai and A2 complex constants. 
The real part of Equation P5)) contains the solution with 



+ A2H^ 



r(2) 



(2) 



2A — exp - — 

„ . / z \ 

2A — exp - — 
vw V 2AJ 



(43) 



physical meaning. The functions Hq and Hq^' repre- 
sent waves propagating towards the positive and nega- 
tive z-directions, respectively. By setting A2 = 0, we 
restrict ourselves to propagation towards the positive z- 
direction. In addition, we chose Ai so that ^{vr) = 1 
and (^) = Q at z = 0, hence 



r(2) 



Al 



Yi (2A^) +iJi (2A^ 



( 2A^) Yi ( 2A-^) - Jl f 2A-^) ¥„ ( 2A^) ' 

(44) 

with J and Y the Bessel functions of the first and second 
kind, respectively. 

Apart from the exact expression of Vr given in Equa- 
tion , we also can obtain an expression in the WKB 
approximation, which is instructive to understand the 
wave properties. Then, by using Equation (j25l) with the 
kink velocity given in Equation (|42p. we get velocity 



Vr(z) « Vr(0) exp ( — - | cxp { i2A- 

V4A/ I Wko 



1 — exp ^ 



(41) and according to Equation (PT|) the wavelength is 



z 

2Ay 

(45) 



2tt { z \ 

A(z) = — vkoexpl^— j . (46) 



We note that Equation p5|) is consistent with the con- 
dition of validity of the WKB approximation, since it 
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coincides with the expression obtained with the asymp- 
totic expansion of for large arguments, i.e., large A, 
is use d in Equation (j43l) (see, e.g.. lAbramowitz fc Stegunl 
119721) . We see that density stratification increases both 
amplitude and wavelength of the kink wave as it propa- 
gates along the loop from the footpoint to the apex. In 
the limit of no longitudinal stratification, i.e., A — > oo. 
Equations (H5|) and P5|) simply reduce to 

Vr(z)^Vr(0)exp( i-^z] , (47) 
27r 

A(^) = A= — z;ko- (48) 

which corresponds to the propagating kink wave in a lon- 
gitudinally homogeneous tube with constant kink veloc- 
ity Dko- Both the amplitude and the wavelength are con- 
stants in this case. 




0.0 0.1 O.Z 0.3 0.4 0.5 

z/L 



Figure 2. Kink wave radial velocity perturbation (solid line), as 
a function of z/L in a longitudinally stratified tube. The different 
lines correspond to l/R = (solid), l/R = 0.1 (dotted), and l/R = 
0.2 (dashed), fn all cases, x = 3, Pio/papox = fO, and ujL/v]^q = 

too. 




Figure 1. Kink wave radial velocity perturbation (solid line), as 
a function of 2/L in a longitudinally stratified tube. In this case, 
the transverse transitional layer is absent so there is no resonant 
damping. Th e symbols correspond to the WKB approximation 
(Equation I I25II ). The dashed lines outline the envelope exp (z/4A). 
Results for X = 3, Pio/Papox = 10, and ujL/v]^q = 100. 

We compare in Figure [1] the exact Vr{z) given by Equa- 
tion with the expression in the WKB approximation 
(Equation (|45l) ) for a particular set of parameters. We 
see that there is a very good agreement between both 
Equations. The effect of longitudinal stratification on 
both the wavelength and the amplitude is clearly seen in 
Figure [TJ 

4.1.2. Case with Resonant Damping 

Here, we include the effect of damping by resonant 
absorption. As an exact solution of Equation ([T^ is very 
difficult to obtain, we use the expressions derived in the 
WKB approximation. First of all, the expression for A(z) 
is the same as given in Equation (|46|) in the case without 
resonant damping. Now we compute the damping length, 
Ld{z), due to resonant absorption using Equation (j33p 
with the kink velocity given in Equation (|^^ . We obtain 



Lu{z) = 27r^^^l«koexp (-1 j . (49) 
mix — la; V2A/ 

We note again that the ratio Lu{z) / X{z) is independent 
of z. 

Finally, using Equations (|37p and ((38|) we obtain the 
expression for Vr{z) and A{z), namely 



Vr{z) w A{z) exp < i2A- 



(50) 



with A{z) the amplitude as a function of z given by 
^(z) = t;,(0)exp(^) 



X exp 



1 m ^ y — 1 . UJ 
A 



1 — exp I 



2AJ 



In Equation ()51|) we see that there are two competing 
effects that determine the amplitude of the propagating 
kink wave. On the one hand, density stratification causes 
the amplitude to increase with height; on the other hand, 
resonant absorption damps the kink wave, so its effect 
is to decrease the amplitude. Whereas the increase of 
the amplitude due to longitudinal stratification is inde- 
pendent of UJ, the damping by resonant absorption does 
depend on uj. Thus, we can anticipate that for a par- 
ticular, critical frequency, Wcrit, the amplitude may be a 
constant independent of z. However, due to the func- 
tional dependence of Equation (|51|) on z, we clearly see 
that the amplitude cannot remain constant for all z. De- 
spite this fact, it is still possible to give an expression for 
the critical frequency for which the amplitude is a ap- 
proximately constant for small z only. To do so, we ap- 
proximate 1 — exp [— z/(2A)] w z/(2A) in Equation ((5T|) . 
Then, we obtain 



2 TO ; X - 1 A ' 



(52) 



For UJ < ojcrit longitudinal stratification is the dominant 
effect and the amplitude increase with z. On the con- 
trary, for u! > Wcrit damping by resonant absorption is 
more efficient and the amplitude decreases in z. We must 
bear in mind that Equation (j52[) is only valid for small 
z, i.e., close to the footpoint of the loop, but gives us 
some qualitative information about the behavior of the 
amplitude depending on the value of the frequency. 

In Figure [5] we plot Vr{z) for fixed frequency and lon- 
gitudinal inhomogeneity scale height, but for different 
values of l/R. In the case l/R = there is no reso- 
nant damping. As l/R increases, the wave amplitude 
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decreases as a result of resonant absorption. However, 
we see that the wavelength is independent of the value 
of l/R. This means that, in the TTTB approximation, 
resonant absorption does not affect the value of the wave- 
length, which is exclusively determined by the longitu- 
dinal inhomogeneity scale height. This is an important 
result from the observational point of view, because the 
wavelength carries information about longitudinal strat- 
ification only (Equation (I^Tj) ). while the amplitude is 
influenced by both longitudinal and transverse density 
profiles (Equation (|38)) ). 

Figure [3] displays the kink wave amplitude computed 
from Equation (j5ip for different combinations of parame- 
ters. In Figure Efa) we determine the effect of l/R when 
the remaining quantities are kept constant, while Fig- 
ure [3Kb) and (c) show the effect of w and Pio/papox, re- 
spectively. Both l/R and ui have a similar effect on the 
amplitude as both of them control the efficiency of the 
resonant damping in front of the increase of the ampli- 
tude due to stratification. When both l/R and u) get 
larger, the effect of resonant damping becomes stronger. 
On the contrary, pio/papex controls the role of longi- 
tudinal stratification through A (Equation (|41I) V As 
Pio/Papcx grows (A decreases), the increase in amplitude 
due to stratification becomes more important. For the 
set of parameters used in Figure EKb) , the approximate 
critical frequency for a constant amplitude according to 
Equation ([5^ is uJcritL/vkQ w 100. We see that, as ex- 
pected, for the critical frequency, amplitude is constant 
for small z, but then increases slightly with larger z. We 
stress again that, for fixed frequency, the amplitude is 



determined by the combined effect of longitudinal strat- 
ification and transverse inhomogeneity, while the wave- 
length is only affected by longitudinal stratification. 

4.2. Time- dependent Numerical Simulations 

Here the aim is to show how the analytical results de- 
rived in the previous sections for the eigenmode problem 
are related to the time-dependent problem. In addition, 
we go beyond the TT and TB approximations used in 
the derivation of the analytical expressions, allowing us 
further insight. To study the problem of linear propagat- 
ing waves. Equations ([l} and ([2]) are numerically solved 
in cylindrical coordinates. Since we are only interested 
in the kink mode we assume an azimuthal dependence of 
the form exp{imip) with m = 1, and the equations are 
integrated in the radial and longitudinal coordinates (r 
and z in our case, respectively). A driver at the foot- 
point of the loop (z = 0), with a single frequency w, is 
introduced to excite propagating waves along the tube. 
The spatial form of the driver in the radial direction is 
chosen, to simplify things, to be the eigenmode of a loop 
with no resonant layer. This enables us to mainly excite 
the resonantly damped modes, while the direct excitation 
of Alfven modes is very weak. Thus, the driver at z = 
is implemented through the periodic variation in time 
and the radial dependence of the kink mode eigenfunc- 
tion s. The simulations are done with the code MoLMHD 
(seeH ona et all 120101 for further details about the nu- 
merical method) in a uniform grid of 1000 x 100 points 
in the r-z plane, and in the range < r/i? < 10 and 
< z/ R < 100. For convenience, in the code time is 
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expressed in units of ta — R/vmo, with UAio the internal 
Alfven velocity at z = 0, so that the dimensionless fre- 
quency is WTA- Note that in the previous sections, the 
dimensionless frequency has been written as WTk, with 
— i/wj^o- The relation between the two time scales is 

TA = T^R/L^2xl{l + x)- 

Hereafter, we concentrate on the analysis of the two 
main variables, Vr and Wi^, i.e., the radial and azimuthal 
velocity perturbations, respectively. Figure U shows the 
two-dimensional spatial distribution of the two velocity 
components, Vr and v^, at a given time and for a partic- 
ular set of parameters. The driver excites propagating 
waves that move upwards along the tube with a charac- 
teristic wavelength and amplitude. The amplitude inside 
the tube decreases with height for both and , while 
the v^p component shows a significant increa se in z within 
the in homogeneous layer (see the work by iPascoe et al.l 
I2OIOI for a longitudinally homogeneous loop). In Fig- 
ure El we perform two cuts along the z-direction of Vr at 
r = and aX r = R. The radial velocity perturbation 
at the center of the loop decreases its amplitude in z due 
to the combined effect of longitudinal stratification and 
resonant absorption. On the contrary, the amplitude of 
Vp at r = R grows in z because of the conversion from 
radial to azimuthal motions that takes place in the in- 
homogeneous layer due to the Alfven resonance. As a 
consequence of this process more energy is concentrated 
in the layer in the form of azimuthal motions as the wave 
propagates upwards along the loop. 

First, we perform some test cases in order to compare 
the full numerical solutions with the analytical expres- 
sions in the TTTB and WKB approximation. We start 
with an example that satisfies both the TTTB and WKB 
approximations. In this example the frequency of the 
driver is wta = 0.36. The transitional layer thickness is 
set to l/R = 0.1, i.e., we consider a a thin layer, and the 
stratification scale height is K/R = 200. Waves with a 
wavelength much longer that the tube radius and much 
shorter than the stratification scale height are excited. 
Therefore, we are in the TTTB and WKB regimes. In 
Figure ini we plot the z-dependence oi Vr at r = 0. The 
radial velocity perturbation shows an almost constant 
amplitude in z, while the wavelength increases slightly 
with height. In this particular example, the effect of lon- 
gitudinal stratification and resonant absorption on the 
amplitude cancel each other. The numerical curve is 
in quite good agreement with the analytical expression 
given by Equation ((50)) . Thus, this simulation shows 
that the assumptions made to derive Equation (|50p are 
well justified when R <^ A ^ A. The almost negligible 
differences between the numerical and analytical results 
become smaller as we take smaller values of R and larger 
values of A. 

Next, we show how the results are modified when we 
depart from the valid regimes of the TTTB and WKB 
approximations. In a second example (see Figure [71) the 
driver frequency is wta = 0.72 and the rest of parame- 
ters are the same as in the previous case, but now the 
excited kink wavelength is almost twice smaller than be- 
fore and the assumption of TT may be compromised. 
The comparison of the numerical result and analytical 
TTTB expression indicates that for large z the analyt- 
ical formula overestimates both amplitude of oscillation 



and wavelength. In our third example, we return to the 
case with wta = 0.36 but we use a larger value of the 
transitional layer thickness, l/R = 0.4. Then, we are 
not in the TB regime. The results in this case (see Fig- 
ure [51) point out that the analytical expressions predict a 
stronger damping of amplitude and a larger increase of 
wavelength with z in comparison with our numerical re- 
sults. In our last case, we test the WKB approximation 
by considering that the wavelength of the excited kink 
mode is not much shorter than the stratification scale 
height. First we take in Figure [9l the same parameters 
as in Figure [6lbut with K/R — 50 and wta — 0.36. Now 
the wavelength and the stratification scale height are of 
the same order, approximately, and we are not in the 
WKB regime strictly. Next in Figure [10] we use exactly 
the same parameters but wta = 0.18, so that the wave- 
length is longer than the stratification scale height and 
the condition of applicability of the WKB approximation 
is violated. However, both Figures [3] and [TUl show that 
the full numerical solutions and the WKB results are in 
remarkable good agreement. 

The numerical simulations performed in this Section 
tell us that the analytical expressions derived in the 
TTTB and WKB approximations are very accurate when 
kink wavelength is in the range i? <C A ^ A and l/R <C 1. 
We also obtain a very good agreement when the condi- 
tion A ^ A is relaxed. For thick layers, i.e., outside the 
TB approximation, the analytical expression for the am- 
plitude gives a slightly stronger damping than in the full 
numerical case, although the wavelength remains approx- 
imately correct. If instead we depart from the TT ap- 
proximation, the analytical approximation of amplitude 
and wavelength gives slightly larger values in comparison 
with the numerical results. However, as can be seen in 
Figures iTl-fTOl the differences between the numerical sim- 
ulations and the analytical predictions are remarkably 
small even in the worst case scenarios studied here. 

5. IMPLICATIONS FOR SOLAR ATMOSPHERIC 
MAGNETOSEISMOLOGY 

The theoretical results of this paper show that deter- 
mining the longitudinal and transverse inhomogeneity 
length scales of solar waveguides from observations of 
propagating kink waves is not a trivial matter. As was 
shown in Section 14.1.21 if '^6 can observationally iden- 
tify the wavelength (or phase speed) as a function of z, 
we know this is independent of the transverse inhomo- 
geneity. Therefore, if we can simultaneously estimate the 
wavelength (or phase speed) and amplitude of the kink 
wave as a function of z we will have enough information 
to quantify the possible contribution of resonant absorp- 
tion on attenuating the observed amplitude. However, 
observationally this presents difficult challenges. In the 
WKB approximation, for a flux tube with density strat- 
ification, constant magnetic field and no resonant layer, 
Equation ([25l) shows that the velocity amplitude is re- 
lated to the phase speed by 

Vr{z) cx yjv\,{z). (53) 

Hence if this relationship does not hold for an observed 
wave, the discrepancy could be caused by damping due 
to the presence of a resonant layer, i.e., transverse in- 
homogeneity. At present a possible avenue to for this 
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Figure 4. Snapshot of the velocity perturbations (a) Vr and (b) v^p at a given time (indicated on top of the panels). The driver is located 
at 2: = and its frequency is ujtx = 0.36. In this plot the width of the transverse transitional layer (see dashed lines) is 1/ R = 0.2, x = 3, 
and A/R = 200. Note the different spatial scale in the r and z axes. 
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Figure 5. Velocity perturbations Vr (solid line) at r = and v,p 
(dashed line) at r = _R as functions of the position along the loop 
at a given time (t/rA = 103.4). These results correspond to the 
simulation shown in Figure |4] 



Figure 7. Simulated radial velocity perturbation (solid line) and 
analytical curve (dotted line) at r = as functions of position along 
the loop at t/rA = 103.4. In this plot l/R = 0.1, A/ R = 200 and 
LOTA = 0.72. 
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Figure 6. Simulated radial velocity perturbation (solid line) and 
analytical curve (dotted line) at r = as functions of position along 
the loop at I/ta = 103.4. In this plot l/R. = 0.1, K/R = 200 and 

LOTA = 0.36. 



Figure 8. Simulated radial velocity perturbation (solid line) and 
analytical curve (dotted line) at r = as functions of position along 
the loop at I/ta = 103.4. In this plot l/R. = 0.4, k/R = 200, and 
LOTA = 0.36. 
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Figure 9. Simulated radial velocity perturbation (solid line) and 
analytical curve (dotted line) at r = as functions of position along 
the loop at t/rA = 70.1. In this plot l/R = 0.1, A/i? = 50, and 
ojTA = 0.36. 




Figure 10. Same as Figure |9] but for ojta = 0.18. 



type of magnetoseismology is by using the high spa- 
tial/temporal resolution Solar Optical Telescope (SOT) 
onboard Hinode. Using the Ca II broadband fiher of 
SOT chromospheric kink waves pr opagating in spicu les 
can be analyzed. As was shown bv iHe et al. I (|2009b[ ). it 
is possible to track both the propagating kink wave am- 
plitudes and p hase speeds from the photosphere up to the 
lower corona. iHe et al~l (|2009aB have also shown time 
snapshots of these propagating waves that nicely show 
the relationship between wavelength and amplitude with 
height. Unfortunately, since there are only a few wave- 
lengths detected up to the visible apex of spicules it is 
difficult to detect a clear trend in the variation of wave- 
length with height. Therefore a more practical approach 
of estimating v\^{z) is by trackin g the phase travel time 
of the kink wave as was d one bvlHe et al. I (l2009bD . 

Regarding CoMP data, iVerth et all H2010f ) attempted 
to estimate the damping length of kink waves propagat- 
ing along coronal loops, assuming the plasma was lon- 
gitudinally homogeneous. To interpret CoMP data, as- 
suming there is longitudinal variation in plasma density, 
as stated previously, we need both information about ve- 
locity amplitude and phase speed to see if Equation (|53|) 
holds. Certainly, Tomczyk & Mcintosh (2009,) showed 
that it was possible to estimate v-kiz) from their time- 
dista nce plots (see Figure 5 in iTomczvk fc Mclntoshl 
[2009f l. However, since the velocity amplitudes were small 



(only a few km s) it will be more of a challenge to es- 
timate the longitudinal variation of this quantity. 



We must also mention the physical limitation of using 
a model fl ux tube with constant magnetic field. It was 
shown by iRuderman et al.l ()2008D that if the magnetic 
field is varying longitudinally (still with density strati- 
fication), in the thin tube approximation Equation (j53p 
becomes 

Vr{z) cx i?(z)Vwk(z), (54) 

where R{z) is th e flux tube r adius . Indeed, Equation 
was used bv IVerth et al.l ()201lD as the basis for in- 
terpreting propagating kink waves in a spicules to deter- 
mine chromospheric flux tube expansion, since both Vr{z) 
and Uk(.z) were estimated from observation. It can be 
seen from Equation (|54|) that if there is a resonant layer 
present, assuming ideal MHD will cause the expansion of 
flux tube with height to be underestimated. Of course. 
Equation (l54t can also be applied to CoMP observations 
of propagating kink waves to estimate the variation of the 
magnetic field along coronal loops since in the thin tube 
flux conservation is given by B{z) oc 1/R^{z). Magneto- 
seismology of post-flare standing kink waves in coronal 
loops, allowing for both longitudinal variation in mag- 
netic field strength and plasma density was first done by 
IVerth et al.l (|2008[ ). therefore it would be interesting to 
see if propagating coronal kink waves could be exploited 
for a similar purpose. 

6. DISCUSSION AND CONCLUSIONS 

In this paper we have investigated the propagation of 
resonantly damped kink MHD waves in both transversely 
and longitudinally inhomogeneous solar waveguides. The 
present work extends the previous investigation by TGV, 
where longitudinal inhomogencity was not taken into ac- 
count. By using the WKB method, we have derived gen- 
eral expressions for amplitude and wavelength of kink 
waves in the TT and TB approximations. Variation of 
wavelength along the magnetic flux tube only depends on 
longitudinal stratiflcation. However, wave amplitude is 
affected by both transverse and longitudinal inhomogenc- 
ity. The kink mode is damped by resonant absorption in 
the Alfven continuum due to transverse inhomogencity 
and the damping length changes along the waveguide 
since it is dependent on the longitudinal density proflle. 
However, it is important to note that the ratio of damp- 
ing length to wavelength along the magnetic flux tube 
is constant. It is obvious that the effect of longitudinal 
inhomogencity can both amplify or attenuate the wave, 
depending on whether the density decreases or increases 
towards the direction of wave propagation, respectively. 
In the case of radial and longitudinal inhomogencity we 
have shown that the variation of kink mode amplitude 
along magnetic tubes is determined by the combined ef- 
fect of both these mechanisms. 

We have performed an application of the theory to the 
case of driven kink waves in stratified coronal loops prop- 
agating upward from the loop footpoints. In this model 
the density decreases exponentially along the flux tube, 
and so resonant absorption and longitudinal stratiflca- 
tion have opposite effects on the kink wave amplitude. 
The efficiency of resonant absorption as a damping mech- 
anism depends on the frequency. For frequencies larger 
than an approximate critical value, the net effect is the 
decrease of the amplitude in z. The opposite effect, 
i.e., increase of the amplitude, takes place for frequencies 
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smaller than the approximate critieal one. These analyti- 
cal results have been checked using full numerical simula- 
tions beyond the TTTB and WKB approximations. The 
analytical expressions and the full numerical solution of 
driven waves are in good agreement even when the re- 
quirements of the TTTB and WKB approximations are 
not strictly fulfilled. 

The analytical expressions for the wavelength and the 
amplitude derived in this paper have direct implica- 
tions for solar magnetoseismology of, e.g., coronal loops, 
spicules, and prominence threads, allowing the equilib- 
rium model to have realistic variation in plasma den- 
sity in both the longitudinal and transverse direction to 
the magnetic field. This is an important step forward 
for exploiting the observations of ubiquitous propagat- 
ing kink waves to probe the plasma fine structure of solar 
atmosphere by implementing magnetoseismological tech- 
niques. 
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